Probing Nonequilibrium Fluctuations through Linear Response 
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Linear response analysis in the nonequilibrium steady state (Gaussian regime) provides two in- 
dependent fluctuation-response relations. One, in the form of the symmetric matrix, manifests the 
departure from the equilibrium formula through the quantity so-called irreversible circulation. The 
other, in the anti-symmetric form, connects the asymmetries in the fluctuation and the response 
function. These formulas represent characteristic features of fluctuations far from equilibrium, which 
have no counterparts in thermal equilibrium. 
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It is now well recognized that the macroscopic response 
property of a system near thermal equilibrium is closely 
related to the dynamical fluctuations of its constituent el- 
ements in the scale down to the mesoscopic to molecular 
size. This striking connection is formulated in terms of 
the linear response theory and known as the fluctuation- 
dissipation theorem (FDT) [l| . The theorem is a conse- 
quence of the microscopic reversibility and directly linked 
to the symmetry of the response function in the system 
with multiple degrees of freedom, i.e., Onsager's recip- 
rocal relation for the transport coefficients Q. These 
concepts play a crucial role for our understanding the 
dynamic hierarchical structure of nature. 

Away from equilibrium, however, the FDT is generally 
no longer valid and associated asymmetries show up as 
a key feature of the nonequilibrium state. This state- 
ment is true even for the system characterized by Gaus- 
sian fluctuations with the broken time-reversal symme- 
try, in which the absence of the detailed balance does 
not allow the characterization of fluctuations from the 
measured response function. Indeed, recent several stud- 
ies have reported the involved fluctuation-response rela- 
tion in nonequilibrium systems with negligible nonlincar- 
ity, in which the equilibrium FDT formula is broken in 
such a way that the correlation and the response func- 
tion evolves differently with time. Examples include the 
driven colloids [3|, dissipative systems, i.e., granular ma- 
terials 0, M, and some phase-ordering systems under 
shear flow |6j. Here the origin of the complexity arises 
not from the nonlinearity but from the coupling between 
different degrees of freedoms. 

The aim of the present paper is to elucidate the 
fluctuation-response relation in nonequilibrium Gaussian 
regime by focusing on asymmetries in cross correlations. 
We shall show that there is no complication due to the 
different temporal dependence in a matrix representation 
and all the nonequilibrium effects appear as the FDT ra- 
tio matrix, which consists of the intensity of the noise 
and the so-called irreversible circulation of fluctuation as 
a manifestation of the violation of the detailed balance. 



To survey the problem under consideration, we first 
employ a simple polymer model under shear flow, for 
which both the correlation and the response function can 
be calculated easily. We then proceed to the general ar- 
gument based on the linear response analysis applied to 
the nonequilibrium steady state, in which the dynam- 
ics of fluctuations obey Gaussian statistics. The char- 
acteristic of nonequilibrium fluctuation-response relation 
is nicely demonstrated by decomposing it into symmet- 
ric and anti-symmetric parts. In particular, the anti- 
symmetric part of the response concerns the deviation 
from the reciprocal relation and one can prove the exact 
relationship between it and the nonequilibrium compo- 
nent of the fluctuation. We then argue that the results 
persist even to the nonlinear dynamics provided that the 
fluctuation around the secular motion is Gaussian, as is 
usually expected for macroscopic systems. 

Polymer in shear flow: Let us consider two beads 
connected by a harmonic spring in the thermal bath of 
the temperature T. Although simple, this model, called 
dumbbell model, enables one to capture basic rheological 
properties of polymer solutions [J. Being placed in, for 
instance, a shear flow, the model provides one of the sim- 
plest examples, in which the detailed balance condition 
is violated. The extension to more realistic model with 
internal modes is straightforward by introducing Rouse 
modes. Each bead bears a dipole, or simply the electric 
charge ±q of opposite signs at both ends so as to be ma- 
nipulated externally by the time dependent electric field 
E{t). The equation of motion for the end-to-end distance 
x = (.t, y, z) of the dumbbell in flow field is given by 



jXj) = -kXi+V3i(t)+&\t) 



(1) 



where 7, k are the friction and spring constant, which 
set the relaxation time To = 7/fc, f^(t) = 2qE(t) 
is the probing force and w(t) is the Gaussian dis- 
tributed random force with zero mean and the variance 
(Wi(t)wj(t')) = 27T%c$0 - f'). The bracket denotes 
the statistical average in the steady state and the re- 
peated indices imply the summation. The Boltzmann 
constant is set to be unity throughout the paper. We 
set the direction of the flow and the velocity gradient 
to be x and y axes, respectively, i.e., the velocity gra- 



* E-mail: sakauc@scphys.kyoto-u.ac.jp 



dient tensor Ki 



7^12 • Since the vorticity (z) di- 
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FIG. 1: Time evolutions of the FDT ratio (x component) 
The ratio is sometimes defined using the sum of all com- 



ponents; ri™ T) (t) = 



asymptotically approach the static values ri^ DT ^ (t) 
W 2 /2) and T s ( ™ T) (i) T(l + W 2 /(2(2 + IV)). 



Both 
T(l + 



rcction is not affected by the flow, we concentrate on 
the dynamics in x — y plane. Equation. ([I]) is easily 
solved, leading to the steady state correlation function 
Cij(t — t') — (xi(i)Xj(t')) in the absence of the external 
field; 



C(t) = je~^ 



ill 

2 



1 



(2) 



where the dimensionless shear rate W — 'jtq is intro- 
duced. The response property is examined by switching 
on the constant electric field at time t — to and following 
the average evolution of the dumbbell toward the new 
stable state; 



(xi{t-t )) p = R ij (t-t )f j 



(p) 



(3) 



where the subscript p indicates the presence of the probe 
force and Rij(t) = J Q Xij(s)ds is the integrated response: 

R(t) = T 1 \ ro ^ J (4) 

« \ 1 - e" / 



At equilibrium (IF 
through FDT; 



0), these quantities are connected 



TX ii {t) = --C ij {t) 



(5) 



or equivalently TRij(t) = CVj(O) — Cij(t). Now, a shear 
flow breaks the detailed balance and the nonequilibrium 

effect is often quantified by the FDT ratio T^ BT '(t) = 
-dX a p/dC a p. In Fig. {TJ, we plot the FDT ratio of the 
dumbbell under shear. The same.or similar plots were 
reported in various systems 0, 0, @] . Individual compo- 
nents of the response and the correlation evolve differ- 
ently with time, thus, the FDT ratio is time-dependent 
in short time and approaches to the static value (> T) 



at longer time t >> To. Our crucial observation is that 
the relation can be cast into the following matrix form: 



(6) 



where the matrix &ij has the dimension of temperature 
and expressed as 



1 —W. 



(7) 



with Ely = [Cij(0)] and the static susceptibility Xij(0) 
(see eq. (fT3"]l). Although compact, this formula is yet 
useless because of the presence of the involved quantity 
Qij . This matrix, while reduces to the bath temperature 
(thus, scalar quantity), i.e., 0^ = T5ij at equilibrium, is 
responsible for all the nonequilibrium effect, which makes 
eq. (J6]) distinct from eq. ([5]). We shall explore the mean- 
ing of this factor in the general linear stochastic system 
discussed below. 

General model: Let us consider a classical system (with 
n gross variables) coupled to a heat bath at tempera- 
ture T. In the Markovian level of description, the dy- 
namical variable x — {x\, X2, ■ • ■ , x n ) obeys the following 
Langevin equation; 



Xi(t) = -KijXjit) + &(t) + v^\t) 



(8) 



where the matrix K^j represents the regression to the 
secular motion, £(t) is a random force and a weak per- 
turbation (t) is switched on when probing the re- 
sponse property of the system. A random force is as- 
sumed to be Gaussian white noise with zero mean and 
the variance (*')) = 2Dij8(t - f). The regres- 

sion matrix can be generally decomposed into two parts 
KiA = + . The first is associated with the con- 
servative force K^'x 



j — Lij VjV(x) where is a mo- 
bility matrix and V(x) = (l/2)Uij X%Xj IS Ch potential en- 
ergy. The matrices Ly and f/y are both symmetric. The 

(2) 

second part , on the other hand, represents the con- 
tribution from non-conservative forces, which drives the 
system out-of-equilibrium. The steady state distribution 
function is 



P(x) 



exp [-(p(x)\ 
J dxexp [—(j){x)] 



(9) 



with a quadratic generalized potential (f>(x) — 

(l/2)Sy XiXj . 

Now we switch on the probe force f- p \t) = LT^Vjty). 
The linear response of the system is described by 



( Xi (t)) p = / dt'X^t -t')&\t') 



(10) 



where Xij (t) is a response function. If one keeps applying 
a constant force f( p ' for an enough time, the system will 
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eventually settle in a new steady state with the average 
and the distribution function given, respectively, by 



tip) 



and 



P p (x) 



x i ) p = X ij (0)f« 

exp [-(f>p(x)] 



J dx exp | 



(11) 



(12) 



where the static susceptibility is defined by (footnote [1]) 



dt X ik (t) = K-/L 3k 



(13) 



By expanding 4> p = (l/2)Y> lj (x i - (xi) p )(xj - (xj) p ) 
and retaining only a linear term in f p , one obtains 

<j> p {x) = 4>(x) - XiMVitWf?* + 0({/ lp) } 2 ) (14) 

Putting eq. (fT4]) into eq. (fT2| gives the following relation 
between P(x) and P p (x) in the linear response regime: 



P p (x) ~ P(x) l + Xy(O)Vi0(f)/ 



(p) 



P(x) 



1 + er 1 ^/. 



ip) 



(15) 



with the matrix Qi k = [Y,ijXj k (0)] 1 introduced in 
eq. ©. 

To examine the dynamic response, we then turn off the 
probe force, say, at t = 0, and follow the relaxation to 
the original steady state. During this relaxation process, 
the average is expressed as 



< X(t) > re lax= J dx xP(x,t) 



(16) 



The probability distribution function at time t is con- 
nected with the initial (t = 0) one through 

W ./rf«W,D (") 

where G(x,x';t) is a propagator of the Fokker-Plank 
equation corresponding to eq. ([5]). By noting P(x, 0) = 
P p (x) and using eq. (fT5j) . one can transform eq. (fTBj) as 
follows 



< x(t) >relax = I dx J dx' X G{x,x']t) P p {x') 

= Qij(t)f^ (18) 

with the relaxation function Qik{t) — f°° ds Xn~(s) = 
Cij(t)Qj£. After differentiation with respect to time, 
one arrives at the fluctuation- response relation: eq. ©. 



[1] Note that the complex admittance is Xik(u)) = (Kij 
iio&ij) -1 Ljk from eq. 1(5} • 



To establish a clear connection between the response 
and nonequilibrium fluctuations, let us analyze the quan- 
tity Qij. The essential feature of the nonequilibrium 
steady state is the presence of non-vanishing probabil- 
ity current. In Gaussian regime, this corresponds to a 
circulating flux, which can be invoked by representing 
the steady state covariance in the form: 



C ik (0)=Kr j 1 (D jk +Sl jk ) 
where the anti-symmetric matrix 



(19) 



(20) 



represents the deviation from the Onsager reciprocity for 
the kinetic coefficient and called irreversible circulation 
of fluctuation 0]. The steady state probability current 
ji(x) = flijV jP(x) is generally a function of a; and diver- 
genceless (Vdi(x) = 0). Using eqs. ©, HU) and (fT9j) . 
the quantity 0y- is decomposed into two constituents; 



©ifc — (Dj k + Qj k ) 



(21) 



Equation (|21[) clearly shows two requisites for the thermal 
equilibrium: 1) a property of noise Dij = TLij, i.e., the 
relation known as FDT of 2nd kind, and 2) no presence 



of the probability current 11; 



0. 



Now we proceed to express eq. |SJ in a frequency 
space. Let us first Fourier-Laplace transform eq. ©: 
Xij(uj)Qj k = Ci k (0) + iuGi k (uj) ,where we have intro- 
duced Xjj(u) — J °° dt Xij(t) e luJt (complex admittance) 
and Gij (u>) = J °° dt Cij (t) e luJt . Using the property of 
the correlation function Cij(t) — Cji(-t) (as a conse- 
quence of the time translational invariance), the power 
spectrum is expressed as Cy (u>) = dt (t) e lu)t = 

Gj^ui) + Gij(uj) (the symbol A*j denotes the complex 
conjugate of Aij.) Then, it follows that 



A^ s) [X'(a;)©] 



(22) 



(23) 



where A.-^fA 



l Aii + Aji ^d Aj s) [A] = Ay - A j{ 

represent the symmetric and anti-symmetric components 
of the matrix Ay, and prime and double-prime denote 
the real and imaginary parts, respectively. Note that the 
imaginary part of the power spectrum is anti-symmetric 
and connected to the breaking of time-reversal symmetry: 



dt [Ci 3 (t) - Cij(-t)] exp (iiot) 



(24) 



From now on, we assume the FDT of 2nd kind, i.e., 
— TL^. While this implies that the microscopic 
degrees of freedom representing the heat bath is equili- 
brated, the FDT of 1st kind is not yet realized in the 
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presence of the non-conservatice force ' . From this 
and eq. (p?Tj) , one obtains 

&ik = TS ik + L^Qjk (25) 

By substituting (f2"5) into eqs. (f2"2")l and (f2"3"|) . these are 
written, respectively, as 

TAg>[X»] - -Ag ) [X"( w )L- 1 ^] (26) 

TA^tX'M] + W C£M - -A J ( ; s) [X'(^)L- 1 fi] (27) 

It is evident from these expressions that the irreversible 
circulation makes the fluctuation-response relation quite 
different from that near equilibrium. That is, near equi- 
librium where = 0, the right-hand side in eq. (|26|) is 
zero (which is referred to as the FDT of 1st kind), while 
both terms of the left-hand size in eq. (|2T|) (C'/j (u>) and 

the anti-symmetric part of X^Aui)) equally vanish, i.e., 
a manifestation of the time reversal symmetry and the 
reciprocal relation. 

Equation (|2T|) can be simplified further by taking no- 
tice of the structure of the anti-symmetric part of the ad- 
mittance. Let us here demonstrate it using the simplest 
case with two gross variables. The real and imaginary 
parts of the admittance are generally connected to each 
other through the Kramers-Kronig relation. For exam- 
ple, in the case of a harmonic oscillator (with one degree 
of freedom) in equilibrium: x = —K\X + £(t), the re- 
lation is given by X'tqui = X", with a relaxation time 
To = . In the present case, one can, indeed, show the 
following similar relation for the anti-symmetric part Q : 

A[ a 2 s) [X»0] = -Tur 12 A^ s) [X»] (28) 

with Qij given by eq. (|2"5")) and t\% — Tl | K ] . Combin- 
ing eq. (f2"B")) with eq. (f2"T|) leads to the following relation 
between the anti-symmetric parts of the admittance and 
the Fourier transform of the cross correlation: 

TTi 2 A% l) [it"(u J )}=C'M (29) 

This formula together with eq. (|2"6"|) , which is the central 
result of the present paper, connects the noncquilibrium 
fluctuation, i.e., the part which is odd under the time re- 
versal (eq. (|24p ). to the asymmetry in the response func- 
tion. 



We emphasize that the formulas ([26]), J28|) and ([29]) 
provide us with the procedure to explore the properties 
of the nonequilibrium fluctuations. If the admittance 
X^(lu) (or X"j(u)) is experimentally accessible, one can 
obtain the regression matrix and the mobility ma- 
trix Lij from its limiting behaviours (see footnote [1]): 
X ik (ui — > 0) = K^j Ljk and X^{u — > oo) = iLijcj^ 1 . 
Then, from eq. (|2"5|) and (f2"9"|) one obtains the irreversible 
circulation and the imaginary part of the power spec- 
trum Cjj(u>), respectively, which are the mostly relevant 
quantities to characterize the systems far from equilib- 
rium. Finally, eq. (|26|) together with the information so 
far obtained provides the real part of the power spectrum 
C'ijfa), thus, the complete knowledge of the correlation 
function. Although our argument is valid only for Gaus- 
sian regime, this does not necessarily exclude the appli- 
cability to the nonlinear dynamics. Even when the gross 
variables obey the nonlinear equation of motion, fluctua- 
tions around the secular motion might be still Gaussian, 
which is usually expected for macroscopic system. In 
such cases, all the derived results can be applied without 
any modifications @. 

In summary, we have constructed the linear response 
formalism around the nonequilibrium steady state in 
Gaussian regime. Main results are (I) the simple equa- 
tion, (|29p . connecting the nonequilibrium fluctuation to 
the departure from the reciprocal relation and (II) the 
equation, (f2"6")) . quantifying the FDT violation through 
the irreversible circulation of fluctuation. These can 
be shown to be consistent with other works on the 
fluctuation-response in nonequilibrium steady state, in 
which the effect of nonlinearity is highlightedflG. KLlL [L2| . 
While restricted, the present discussion in Gaussian 
regime enable us to establish the useful formulas elucidat- 
ing connections between essential quantities characteriz- 
ing nonequilibrium steady state, through which one can 
extract the role of the coupling between different degrees 
of freedom away from equilibrium. 
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